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ABSTRACT 

We use the complete, X-ray flux-limited ROSAT Bright Survey (RBS) to measure 
the space density (p) of magnetic cataclysmic variables (mCVs). The survey pro- 
vides complete optical identification of all sources with count rate > 0.2 s" 1 (cor- 
responding to Fx > 2 x 10~ 12 ergcm^s" 1 ) over half the sky > 30°), and de- 
tected 6 intermediate polars (IPs) and 24 polars. If we assume that the 30 mCVs 
included in the RBS are representative of the intrinsic population, the space den- 
sity of mCVs is 8^2 x 10~ 7 pc~ 3 . Considering polars and IPs separately, we find 
Ppoiar = 5^2 x 10~ 7 pc~ 3 and pip = 3t 2 x 10~ 7 pc~ 3 . Allowing for a 50% high-state 
duty cycle for polars (and assuming that these systems are below the RBS detec- 
tion limit during their low states) doubles our estimate of p po iar and brings the total 
space density of mCVs to 1.3~t ' A x 10~ 6 pc" 3 . We also place upper limits on the sizes 
of faint (but persistent) mCV populations that might have escaped detection in the 
RBS. Although the large uncertainties in the p estimates prevent us from drawing 
strong conclusions, we discuss the implications of our results for the evolutionary re- 
lationship between IPs and polars, the fraction of CVs with strongly magnetic white 
dwarfs (WDs), and for the contribution of mCVs to Galactic populations of hard X- 
ray sources at Lx > 10 31 ergs _1 . Our space density estimates are consistent with the 
very simple model where long-period IPs evolve into polars and account for the whole 
short-period polar population. We find that the fraction of WDs that are strongly 
magnetic is not significantly higher for CV primaries than for isolated WDs. Finally, 
the space density of IPs is sufficiently high to explain the bright, hard X-ray source 
population in the Galactic Centre. 

Key words: binaries - stars: cataclysmic variables, - X-rays: binaries - methods: 
observational, statistical. 



1 INTRODUCTION 

Cataclysmic variables (CVs) are interacting binary stars in 
which a white dwarf (WD) accretes matter from a low-mass, 
Roche- lobe-filling companion. In around 20% of known CVs, 
the magnetic field of the WD is sufficiently strong to control 
at least the inner part of the accretion flow; these magnetic 
CVs (mCVs) are divided into two classes, namely polars and 
intermediate polars (IPs). The defining property of a polar 
is that the WD is locked in synchronous rotation with the 
binary orbit, while IPs have WD spin periods typically much 
shorter than the orbital period (P or b)- 

In many ways, the formation and evolution of mag- 
netic and non-magnetic CVs is thought to be similar. Both 
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types of systems form through common envelope (CE) evo- 
lution, evolve first from long to short P orb because of an- 
gular momentum loss (AML), and eventually experience 
period bounce, when the thermal time-scale of the donor 
becomes longer than its mass loss time-scale. In fact, the 
main proposed difference between the evolution of mCVs 
and non-magnetic CVs affects only the polars. In these sys- 
tems, magnetic braking (MB), which is thought to be the 
dominant AML mechanism for most CVs above the period 
gap, is likely to be suppressed ( e.g. iLi fc Wickramasinghd 
1998; iTownslev fc Gansickd 12009). The Porb distribution of 
mCVs is broadly in line with these ideas: if polars and IPs 
are considered jointly, their P oro distribution is very similar 
to that of non-magnetic CVs, showing both a period gap in 
the range 2 hr < P or b ^ 3 hr and a period minimum at around 
P orb ~ 80min (see Fig. [TJ. 

However, there are at least three important open ques- 
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Figure 1. The orbital period distribution of all CVs (fine his- 
togram), and mCVs (bold histogram). Cumulative distributions 
are also shown for polars (dashed histogram) and IPs (dotted 
histogram). Almost all IPs are found at long P or h, while most 
polars are short-P or (, systems. Although the period distribution 
of mCVs shows a period gap, it is less pronounced than for non- 
magnetic CVs. The periods used here are taken from version 7.18 
of the catalogue of Ritter & Kolb (2003). 



tions concerning the formation, evolution and Galactic 
abundance of mCVs: 

(i) Is there an evolutionary relationship between IPs and 
polars? 

(ii) What is the intrinsic fraction of mCVs amongst the 
general CV population, and can this be reconciled with the 
incidence of magnetic WDs in the isolated WD population? 
(iii) Do mCVs dominate the total Galactic X-ray source 
populations above Lx ~ 10 31 ergs -1 ? 

In principle, all of these questions can be addressed empir- 
ically, but this requires reliable measurements of the space 
density of IPs and polars. The goal of the present work is to 
provide such measurements. 

The plan of this paper is as follows. In Section [5] we will 
provide some additional background on the questions listed 
above, and discuss previous attempts to measure the space 
densities of polars and IPs. Our flux-limited, X-ray-selected 
mCV sample is presented in Section [3] with distance and 
X-ray luminosity estimates discussed in Section U and [S] In 
Section [(J we describe how we use this sample to calculate 
space densities, and present our results. The results are dis- 
cussed in Section [3 and, finally, we list our conclusions in 
Section [3 



2 CONTEXT 

2.1 The relationship between intermediate polars 
and polars 

It has been known for a long time that most IPs are found 
above the period gap and most polars below (see Fig. [T]). 
This immediately suggests that IPs may evolve into po - 
lars (|Chanmugam fc Ravlll984l ; lKing. Frank fc Rittedl 1985t h 
This is a physically appealing idea. Whether the WD in a 
given system will synchronize depends on the strength of 
the WD's magnetic field, the orbital separation, and the 
mass-transfer rate (M) from the secondary star (with large 



magnetic field, small orbital separation, and low M favour- 
ing synchronization). Since MB drives much higher mass- 
transfer rates above the period gap than gravitational radi- 
ation (GR) does below, it is plausible that many accreting 
magnetic WDs may only achieve synchronization once they 
have crossed the period gap. 

There is, however, a serious problem with this attrac- 
tive scenario. If most IPs evolve into polars, and most 
polars are the descendants of IPs, their field strengths 
should be comparable. Yet, empirically, the magnetic 
fields of the WDs in IPs (B IP < 10MG) are sys- 
tematically weaker than those of the WDs in polars 



(B volar ~ 10 - 100MG) (e.g. 


Chanmugam fc Ravi 1984 ; 


Wickramasinghe, Wu & Ferrarid 


1991 ; Lamb & Melialll987 ; 


Schmidt ct al. 19961. l200ll; Norton, Wvnn & Somcrscalcs 


2004 


; SchwoDe. Schreiber & Szkodv 20061; Butters et al. 


2009 


)• 



There are several possible resolutions to this problem. 
Perhaps the simplest (in terms of binary evolution) is that 
the high accretion rates in IPs may partially "bury" the 
WD magnetic fields, so that the observationally inferred 
field strengths for t hese systems a re systema t ically biased 
low (|Cummin 3 120021 '). Alternatively. IPattersonI |l994T ) points 
out that it may not be necessary for all (or even most) long- 
period IPs to evolve into polars in order to account for the 
entire short-period polar population, since the evolutionary 
time-scales are so much shorter above the period gap than 
below. This idea still requires an explanation for the fate of 
the remaining long-period IPs, although perhaps they sim- 
ply become unobservable as their accretion flows become 
clumpy and their luminosity is pushed into the EUV band 
|Wickramasinghe. Wu fc Ferrariol ll991). Finally, the oppo- 
site view would be that the short-period polar population 
may be dominated by systems born below the period gap. 
But in this case, there would have to be a direct link be- 
tween the WD field strength and the orbital perio d after 
emergence from the CE phase (e.g. iTout et al.l l2008). 

One way to shed light on the relationship between IPs 
and polars is via their respective space densities. For ex- 
ample, if all long-period IPs evolve into short-period po- 
lars, and all short-period polars are the progeny of long- 
period IPs, then their space densities should be propor- 
tional to the evolutionary time-scale associated with these 
phases. In this particular example, we would predict that 

Ppolar/piP — TGr/tmb » 1. 



2.2 The intrinsic fraction of magnetic CVs 

Magneti c systems make up ~ 20% of the known CV pop- 
ulation jRitter fc Kolbl f2003). At first sight, this is a sur- 
prisingly high fraction, given that the strong magnetic fields 
characteristic of IPs and pol ars ( B > 10 6 G) are found in only 
~ 10 % of isolated WDs (e.g. lKawka et al.ll2007l ; lKiilebi et all 
2009). If these numbers are representative of the intrinsic 
incidence of magnetism amongst CVs and single WDs, the 
difference between them would have significant implications: 
either strong magnetic fields would have to favour the pro- 
duction of CVs, or some aspect(s) of pre-CV evolution would 
have to favour the production of strong magnetic fields. 

A model a long t he latter lines was recently proposed 
by ITout et~aT1 (120081 ). They argue that all strongly mag- 
netic WDs are the product of binary evolution, and that 
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the observed strong fields develop during the CE phase. In 
this picture, mCVs are associated with CE events that pro- 
duce very close binary systems (which are then brought into 
contact by MB or GR), while isolated magnetic WDs are 
produced when two degenerate cores undergo a complete 
merger during the CE phase. 

However, it is actually by no means clear yet that mag- 
netism is really more common in CV primaries than in iso- 
lated WDs. The main problem is that the observed fraction 
of magnetic systems amongst known CVs is almost certainly 
affected by serious selection biases. For example, since mCVs 
are known to be relatively X-ray bright, they are likely to be 
over-represented in X-ray-selected samples. Conversely, po- 
lars, in particular, are relatively faint in optical light (since 
they do not contain optically bright accretion disks), so they 
are likely to be under-represented in optically-selected sam- 
ples. Given that the overall CV sample is a highly hetero- 
geneous mixture of X-ray-, optical- and variability-selected 
sub-samples (which also usually lack clear flux limits), it is 
very difficult to know how the observed fraction of mCVs 
relates to the intrinsic fraction of magnetic WDs in CVs. 



2.3 Galactic X-ray Source Populations 

There have been many attempts to determine the make- 
up and luminosity function of Galactic X-ray source pop- 
ulations in a variety of enviro nments, rangin g from the 



Milky Way as a wh ole (e.g. Sazonov ct al. 2006 



)6), 
39). 



the 
the 



Galactic Centre (e. g. iMuno et all 12004 120061 [200! 
Gala c tic Ridge (e .g. Revnivtsev et al.ll2006l ; iKrivonos et al.l 
20071: iHone et all |2012T ).~ and even globular clusters (e.g. 
Heinke et all I2005T ). Remarkably, in all of these environ- 



ments, mCVs have been proposed as the dominant popu- 
lation of X-ray sources above Lx ^ 10 31 ergs _1 . 

However, in most of the above studies, the breakdown 
of the observed X-ray source samples into distinct popu- 
lations is subject to considerable uncertainty. In general, 
only a small fraction of the sources have optical counterparts 
and/or properties that would permit a clear classification. 
Thus, for the most part, identifications of observed sources 
with physical populations have to rely on gross properties 
(such as X-ray colours and luminosities) and statistical com- 
parisons of observed and expected number counts. The local 
space densities of the relevant populations are arguably the 
most important ingredient in these comparisons. In effect, 
the question being asked is whether the extrapolation of the 
local space density to the environment being investigated 
can account for the observed number of sources seen there. 
In the case of mCVs, such extrapolations are difficult, pri- 
marily because the local space densities are rather poorly 
constrained. 



2.4 The Space Density of Magnetic CVs: 
Previous Work 



2.4.1 Polars 



Estimates for the space d ensity of polars in clude ~ 3.5 x 
10~ 7 pc- 3 dPattersonlll984). ~ 3 x 1 0~ 7 pc~ 3 (IWarnerlll995h . 
~ 5x 10 7 pc~ 3 (ICropperl Il99d ). and 1.3 x lO^pc -3 
jArauio-Betancor et al. 2005). As emphasized particularly 
in the last of these studies, all of these numbers are, strictly 



speaking, lower limits, because they are effectively based on 
the number of systems observed within a particular distance 

CUt-Off; 

iThomas fc Beuermannl (1 19981 ) derive an estimate of 
p po iar — 6.1 x 10~ 7 pc ~ 3 from the sample of polars d etected 
by ROSAT (see also H cucrmann fc SchwopdlT994h . How- 
ever, they also note that, allowing for incompleteness, the 
space density could be as high as 1.2 x 10 -6 pc -3 for short- 
period polars and 1.7 x 10~ 7 pc -3 for long-period polars. 
I Ramsay et al. (2004b) discuss the issue of incompleteness 
in observed samples, particularly that arising from the ex- 
tended low-states exhibited by polars. 



2.4-2 Intermediate Polars 



Fewe r estima tes are available for the space density of 

IPs. IWarnerl (Il995h gi ves p IP ~ 3 x 10" 8 pc" 3 , while 
iRevnivtsev et alj (|2008f ) find p IP ~ 1.2 x 10" 7 pc" 3 , from 
a sample of CVs detected by INTEGRAL. The scarcity of 
space density estimates for IPs is particularly problematic, 
because it is IPs, rather than polars, which are often as- 
sumed to be the dominant bright hard X-ray source popu- 
lation in many Galactic environments (see Section \7.'A\ . 



3 THE FLUX-LIMITED CV SAMPLE 

Given that we do not yet have a useful volume-limited CV 
sample, the most suitable sample to use for statistical stu dies 
is a purely flux-limited one (e.g. iPretorius et al.ll 2007a 



In two previ o us p apers (jPretorius et al.l 120071: 
IPretorius fc Kniggel I2012T ) we presented constraints on 
the space density and X-ray luminosity function of non- 
magnetic CVs, base d on the ROSAT Bright Survey (RBS; 
ISchwope et al.l 12000'). as well as on the deeper, but smaller 
area ROSAT North Ecliptic Po le (NEP) survey (e.g. 
iGioia et al.ll2003l ; iHenrv et alj|2006h . Here, we use the RBS 
to derive the corresponding observational constraints on the 
mCV population. 

While several non-magnetic CVs are known to have very 
faint X-ray luminosities, mCVs are expected to be intrinsi- 
cally brighter in X-rays, because of the almost radial accre- 
tion flow (polars in particular also have higher soft X-ray 
to optical flux ratios than non-magnetic CVs, explaining 
w hy so many where first disco vered by ROSAT; see fig. 2 
of lBeuermann fc Thomas] 1 19931 ). This means that the bright 
flux limit of the RBS is not as important a concern when 
studying mCVs. 

The RBS is a flux-limited part of t he ROSAT All-Sk y 
Survey (RASS; see lVoges et ai1ll999l and lVoges et alJl200Ch , 
consisting of bright (count rate > 0.2 s _1 ), high Galactic 
lati tude (\b\ > 30°) sources. It has complete optical follow- 
up l|Schwope et al. 2000; ISchwope et alj|2002h . and includes 
30 mCVs — 6 IPs and 24 polars. Assuming a 30 keV ther- 
mal bremsstrahlung spectrum, and a 30 eV blackbody spec- 
trum (see Section [5]), both with Nh = 10 20 cm -2 , the 
limiting count rate corresponds to flux limits of Fx > 
3 x 10~ 12 erg cm~ 2 s~ 1 and 2 x 10~ 12 ergcm~ 2 s _1 , respec- 
tively, in the 0.12-2.48 keV band . 

In Table [T] we list the 30 mCVs that make up our 
flux-limited sample. It has been suggested that two more 
CVs detected in the RBS are magnetic; these are TW 
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Pic fsee iMouchet et all Il99ll : iPatterson fe Moulded 1 19931 
iNorton et al.ll2000l ) and V405 Peg l|Thorstensen et aljboogh . 
We do not include these two systems here, since there is no 
conclusive evidence of a magnetic nature for either of them. 
On the other hand, CC Scl was until recently classified as 
a n ormal SU UMa sta r, but has now been shown to be an 
IP (jWoudt et alj|2012h : we therefore include it in our mCV 
sampl^E 



4 DISTANCE ESTIMATES 

We require distance estimates for all mCVs in the sample. 
Where possible, we use published distances, based on either 
trigonometric parallax, or photometric parallax of the WD 
or donor star. In most cases, however, the best distance es- 
timates we can obtain are quite uncertain. This will be the 
most important limit on the precision of our space density 
measurement. The estimates we use are listed in Table [T] 
and we provide more information below. 



4.1 Systems with reliable distance estimates 

4-1.1 Trigonometric 



The IPs EX Hya and TV Col, and the polar EF Eri have 
very reliable distance estimates from parallax measurements 
ilMcArthur et al.ll200ll ; iBeuermann et al.ll2003l ; iThorstenser] 
l2003ft . 



Il993h , we can use the p redicted abso l ute m agnitude of 
the star l|Kniggei 120061 ; iKnigge et ail l201ll ) to find a 
distance of 290±«n° pc, in line with the ~ 250 pc re- 
ported by ISchwope. Thomas, fc Beuermannl (U993). The 
smaller distance estimat e s of Ciardi et al.l Q1998) and 
ISproats. Howell, fc Masonl (1 19961) result from assuming a 
later spectral type. Gansicken 19991 ) also estimated ~ 180 pc 
from low-state UV data, but he notes that there was still 
some accreti on during the observa tion. 

Finally, iReinsch et ail (|l994j)_ estimate a distance of 
300 ± 60 pc for UW Pic, and ISchwarz et all (|l998h find 
620 ± 100 pc for AI Tri, both also from the photometric par- 
allax of the donor. 



4.2 Systems for which only less reliable estimates 
are possible 

For the remaining 22 systems, we have to turn to a less direct 
method of estimating distances. The donor star is expected 
to be responsible for a large fraction of the near- and mid-IR 
flux of a CV. Therefore, even in cases where the contribu- 
tion of the donor star to a system's IR light is not known, 
it is possible to use IR photometry to obtain rough distanc e 
estimates. This is similar to the method of iBailevl (| 198 ll ). 
althoug h we wil l use the p redicted abso l ute do nor magni- 
tudes of lKniggel i|2006h and lKnigge et al.l l|201ll ). Below, we 
will first estimate the average contribution of the donor to 
the near- and mid-IR flux of mCVs, and then use this to 
estimate distances for most of the systems in our samplfl 



4-1.2 Photometric parallax 

Several more systems have reliable distance estimates from 
photometric parallax of one of the stellar components. These 
are based on, e.g., the detection of donor star features in the 
optical or near-IR spectrum (from which the donor's contri- 
bution to the total flux in a given band can be estimated), or 
FUV observations in a low state, where the WD dominates 
the flux. 

Bas ed on donor star features in th e near- IR spectrum of 
DO Dra. lMateo. Szkodv. fc GarnavTc hi (|l99ll ) find a distance 
of 155 ±35 pc. This is very close agre ement with an estim ate 
fro m the FUV detection of th e WD jHoard et al.ll2005l ). 

lArauio-Betancor et al.l (|2005l ) find a distance of 
163^26 for BL Hyi, from a FUV detection of the 
WD (several more estimates, all consistent with this, 
have b een derived from detections of the do nor star ; 
see IVisvanathan. Bessell. fc Wickramasinghel 1 19841. 



Beuermann et aiT ~ 19851 . Glenn et al.l 19941 . Beuermannl 
20001 ). 

HU Aqr is a deeply eclipsing polar; I Ciardi et al.l (|l998l ) 
find K = 15.3(1) in eclipse, and note that the eclipse is prob- 
ably total. Knowing the o rbital period and the spectral type 
of the secondary (M4.5; ISchwope. Thomas, fc Beuermannl 



Although we previously incorrectly included CC Scl in the RBS 
non-magnetic CV sample, it made only a very small contribution 
to our earlier estimate o f the non-magnetic CV space density 
llPretorius fc Knigge 2012), so that the result presented there was 
not significantly affected. It is however important to include CC 
Scl here, because the total space density of magnetic systems is 
smaller. 



4-2.1 Finding the average mCV donor contribution to 
near- and mid-IR flux 

iKniggel (|2006l ) and IKnigge et ail (|201ll ) give the average 
differences between the absolute JHK magnitudes of CV 
donors predicted by their semi-empirical sequence and the 
measured absolute magnitudes of a sample of CVs (consist- 
ing mostly of non-magnetic systems, but including also a few 
mCVs) with parallax distance estimates. This gives an esti- 
mate of the typical donor contribution to the flux in these 
bands, and, together with an apparent IR magnitude, can 
be used to obtain a rough distance estimate (as well as a 
robust lower limit on the distance). 

In order to repeat this for mCVs alone, we use a sample 
of 23 magnetic systems with orbital periods below 6 hours 
and reliable distance estimates. The distances are based on 
trigonometric parallax (for 4 IPs and 5 polars), photometric 
parallax of the WD (5 polars), or photometric parallax the 
donor (1 IP and 8 polars). Not surprisingly, we find that 
the typical donor contribution to the IR light is similar for 
IPs and non-magnetic CVs (systems with discs), and smaller 
for polars. Therefore, we group IPs and non-magnetic CVs 
together, and consider the sample of polars separately. Fur- 
thermore, we decided to omit the polars with published dis- 
tances from photometric parallax of the donoo this means 



2 lAk et ah! ||2008|) present another method of estimating distances 
from near-IR, photometry, and also give distances for several of 
the mCVs in this sample. 

3 Because it is possible that in these systems, the contribution to 
the IR light of, e.g., the accretion flow, is systematically low, given 
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Table 1. The 30 mCVs detected in the RBS, together with their subtype (polar or IP), orbital periods, Galactic latitudes, ROSAT 
PSPC count rates, distances, and X-ray luminosities. The fraction of the total mCV space density contributed by each system 
(Pj/Po) i s given in the 8th column (this ignores all errors, although uncertainties are correctly accounted for later; see Section l6.1jl . 
References are for the classification as a magnetic CV, P or h, and published distances, where available. Note that the distance and 
Lx estimates for CD Ind, IW Eri, CV Hyi, and FH UMa are very uncertain (see Section 14.2.31 and 16. 3. U . 



System 


RBS# 


Type 


P or j,/h 


6 


count rate/s 


1 d/pc 


log(Lx/ergs 1 ) 


Pi 1 Po 


References 


CC Scl 


1969 


IP 


1.402 


—68.7° 


0.28(6) 


nnf . 4-110 

200-70 


31.3(4) 


0.089 


1,2,3 


EX Hya 


1173 


IP 


1.638 


+33.6° 


5.4(3) 


64.5 ± 1.2 


31.59(4) 


0.043 


4,5 


AO Psc 


1914 


IP 


3.591 


—53.3° 


0.37(4) 


oon+180 

,mu_ 120 


32.2(4) 


0.048 


6.7 


DO Dra 


1022 


IP 


3.969 


+44.5° 


0.59(3) 


155 ± 35 


31.4(2) 


0.119 


8,9 


TV Col 


655 


IP 


5.486 


—30.6° 


0.36(4) 


37 °ii5 


rv 00 1 U . U 

32.22^0.09 


0.046 


10,11 


hjl UMa 


713 


IP 


6.434 


+37.4° 


0.57(4) 


' ou _200 


no 9 + 0.1 

13 'J- -0.3 


0.021 


12,13 


CV Hyi 


213 


AM 


1.297 


—50.7° 


0.28(4) 


Kcri+450 
OOU -250 


32.3(6) 


0.013 


14 


T T A TO O O 

V4738 bgr 


1678 


AM 


1.300 


—33.9° 


0.33(4) 


250 -70 


qi O+0.4 

3L8 -0.5 


0.035 


14 


EV UMa 


1219 


AM 


1.328 


+63.1° 


1.86(6) 


cnn+270 
690 -190 


no ( \ 


0.003 


15,16 


GG Leo 


842 


AM 


1.331 


+49.0° 


1.1(1) 


170 -60 


31.5(4) 


0.058 


17,18 


T71T T T TTV T 

bti UMa 


904 


AM 


1.336: 


+48.4° 


0.26(2) 


cnn+480 

590I 270 


00 c\ ( n\ 

32.0(6) 


0.023 


19 


EF Eri 


398 


AM 


1.350 


—57.4° 


6.2(3) 


163 -50 


32 -°±0.4 


0.023 


20,21 


IW Eri 


541 


AM 


1.452: 


-40.6° 


0.23(3) 


270+^30 


31.6(6) 


0.048 


22 


T — 1 T T T TTV T 

EU UMa 


1039 


AM 


1.502 


+76.3° 


3.3(2) 


240 -60 


32.4(3) 


0.010 


18,23,24 


EQ Cet 


206 


AM 


1.547 


—80.9° 


0.34(3) 


270 -90 


31.4(4) 


0.077 


22,25 


V393 Pav 


1664 


AM 


1.647 


—31.3° 


0.90(7) 


o , n +170 
340 -110 


32.4(5) 


0.011 


26 


EG Lyn 


696 


AM 


1.656 


+34.5° 


0.25(3) 


._„+230 
4701 !60 


32.4(4) 


0.010 


22,27 


RS Cae 


599 


AM 


1.699: 


-39.1° 


1.20(8) 


880^30 


33.0(3) 


0.004 


28 


CD Ind 


1735 


AM 


1.848 


-41.4° 


0.38(4) 


350 i 

— 100 


31.7(4) 


0.044 


29,30 


BL Hyi 


232 


AM 


1.894 


-48.6° 


2.8(2) 


163+i* 


31.7(2) 


0.037 


31,32,33,34,35 


EK UMa 


911 


AM 


1.909 


+55.2° 


1.11(5) 


590+ 340 


32.6(4) 


0.008 


36,37,38 


AN UMa 


938 


AM 


1.914 


+62.1° 


1.77(7) 


300+ 150 


31.9(4) 


0.029 


39,40,41 


V1007 Her 


1646 


AM 


1.999 


+33.3° 


0.24(2) 


560±f 8 ° 


32.3(5) 


0.011 


42 


HU Aqr 


1724 


AM 


2.084 


-32.6° 


0.81(7) 


290+n° 


09 9+0.4 

ciZ - Z -0.5 


0.015 


43,44,45 


UW Pic 


658 


AM 


2.223 


-32.7° 


0.73(9) 


300 ± 60 


31.9(3) 


0.028 


46,47 


RX J0859 


734 


AM 


2.397 


+30.9° 


0.23(3) 


5301?™ 


32.6(4) 


0.008 


48,49,50 


CW Hyi 


324 


AM 


3.030 


-45.9° 


0.26(4) 


500±»g 


32.4(4) 


0.035 


22 


RX J1610 


1563 


AM 


3.176 


+37.3° 


0.36(4) 


380±gg 


32.1(5) 


0.050 


22,51 


1RXS J231603 1973 


AM 


3.491 


-58.6° 


1.10(9) 


460+ 170 


32.8(3) 


0.026 


51,52 


AI Tri 


274 


AM 


4.602 


-30.3° 


0.54(5) 


620 ± 100 


32.5(3) 


0.031 


53 
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that our average donor contributions for polars are based on 
a sample of only 10 systems. 

For the two samples of CVs with reliable distance es- 
timates (IPs , together with the non-magnetic CVs used by 
and the sample of 10 polars) separately, we cal- 
culate the average of the difference between the measured 
absolute magnitude of the CVs and the absolute magni- 
tude p redicted for the donor star at that period by [K nigge 
(|2006h . We do this for the 2MASS (|Skrutskie et all 120061 ) 
JHK bands, as we ll as the 2 shortest wavelength WISE 
|Wright et alj l201Ch bands, Wl and W2. On average, the 
donor stars in polars emit, e.g., around 60% of the if -band 
light (compared to 31 %, for the sam ple of of mostly non- 
magnetic CVs used by iKniggel 120061 ) . These results will be 
presented in more detail in a future paper. 

Clearly, a single number cannot describe the donor con- 
tribution to the flux in some waveb and for a l l syste ms and at 
all P rb- However, the sequence of Knigge (2006), together 
with the typical fractions of the near- and mid-IR flux con- 
tributed by the donor, should allow us to predict absolute 
magnitudes to within the errors implied by the CV sam- 
ples used to find these fractions. We take the error in the 
predicted absolute near- and mid-IR magnitudes to be the 
standard deviation of the absolute magnitudes of these two 
samples from the sequence, after it is offset to include the 
typical flux contribution from sources other than the donor. 

A final caveat is that the donor sequence may not be 
applicable to long-period polars, because it is based on the 
masses and radii of donor stars in non-magnetic CVs (which 
have MB above the orbital period gap, while polars perhaps 
never experience MB). We are probably safe in ignoring this 
difficulty, since we use the method for only 3 long-period 
polars, and since no single distance estimate strongly affects 
our space density estimat^S 

4.2.2 Distance estimates relying on the estimated typical 
donor IR flux contribution 

We use the method outlined above to find distance esti- 
mates for 18 of the mCVs in the RBS sample. In doing this, 
we assume that none of these systems are period bouncers, 
or have evolved donor stars (except for EI UMa, where we 
adopt an estimate from the literature that allows for this; see 
below). We include the effect of interstellar extinction in the 
distance estimates when using 2MASS photometry, but ne- 
glect it when using WISE mid-IR photometry (extinction es- 
timates for our sources are discussed in Section[5]below). We 
then find the probability distribution function for the dis- 
tance to each source, assuming Gaussian errors in predicted 
absolute magnitudes, apparent magnitudes, and extinction. 
The distance estimates listed in Table Q] are the median, to- 
gether with the 1-cr confidence interval corresponding to the 
16th and 84th percentile points. We briefly discuss individ- 
ual systems below. 

that spectral features of the donor stars are relatively prominent 
in the optical spectra. 

4 Furthermore, the high mass loss rates of non-magnetic CVs 
above the peri od gap cause the the donor radius to expand by 
at most ~30% l lPatterson et al.ll2005l: lKniggell20u6h . The implied 
~30% distance error is too small to dominate the uncertainty in 
our estimates. 



CC Scl has DN outbursts (as several IPs do), but since 
it is certainly not a typical DN (besides being an IP, its 
superoutbursts have unusually low amplitude and short du- 
ration), we prefer not to base our distance estimate on the 
outburst maximum here. We obtain a distance of 200^].q pc 
from WISE photo metry, slightly sm aller than the estimate 
of 360± 130 pc that lPattersoil l|201lf) finds from the outburst 
maximum. 

The distance of 330±i|g pc that we find for AO Psc 
using WISE photo metry is consisten t with the estimate of ~ 
250 pc obtained bv lPattersonl (|l984h from a relation between 
the strength of emission lines and the absolute magnitude 
of the disc. 

F or EI UMa, we adopt the estimate of iReimer et al.l 
(2008), which is based on the same method we use here. 
Note that it is a tentative estimate because this system has 
P or b in the range where the CV population is expected to be 
dominated by systems with evolved donors (some attem pt 
was made to reflect this in the error; [Reimer et al.ll2008T l. 

EV UMa was in a high state when it w as observed by 
2MAS S; we therefore use K = 1 9.3(2) from lOsborne et all 
(|l994h to obtain d = 690±^ pc. lOsborne et all (|l994h giro 
a lo wer limit of 70 5 pc on the distance to this system. 

iRamsav et al.l (|2004al ) find distances of 50-70 pc and 
60-80 pc for GG Leo and EU UMa, respectively, from UV 
data. They note, however, that these are likely underesti- 
mates, since some of the UV flux might be from sources other 
than the WDs. We find larger distances in both these cases 
(using WISE photometry for GG Leo, and the 2MASS H 
measurement for EU UMa). Our estimate of d = 170^60 P c 
for GG Leo is consiste nt with limits of d > 100 pc (based on 
the optical spectrum; iBurwitz et al.1 Il99ct) and d < 300 pc 
(from the requirem ent that the implied acc retion rate is not 
unreasonably high; iBrinkworth et al.|[2007T l. 

We use WISE data of V393 Pav to estimate d = 
3401! In pc, in agreem ent with the estimate of d ~ 350 pc by 
iThomas et ail (| 19961 ) . 

V4738 Sgr and RS Cae are not detected by 2MASS or 
WISE; we find H = 17.61(9) and 18.1(1), for these two sys- 
tems, from ima ges obtained with the In frar ed Survey Facility 
(IRS F; see e.g. M 

ass fc N agata 2 000l and iNagavam a et all 
2003) at the South African Astronomical Observatory. This 
gives distance estimates of 250^7° pc and 88OI240 P c f° r 
V4738 Sgr and RS Cae, respectively. Note that although 
the r adial velocity curve of RS Cae is aliased (|Burwitz et al.l 
Il996h the uncertainty in period does not significantly affect 
the distance estimate. From the absence of M star features 
in th e optical spectra, lo wer limits of d > 440 pc for RS 
Cae (IBurwitz et al.l [l996T ) and d > 190 pc for V4738 Sgr 
l|Burwitz et al. 1997) hav e been derived. 

Clayton fc Osbornd l| 19941 ) report K = 16.4(7) for EK 
UMa (which is also not detected by 2MASS or WISE), and 
give a lower limit of d > 410 pc. We use their K-h&nA mag- 
nitu de to estimate d = 590^210 P c f° r this system. 

Patterson! (|l984l ) estimates d ~ 400 pc for AN UMa 
(Warner 1995 lists a limit of > 270 pc; see also lLiebert et all 
1982 for some discussion on lower distance limits). We find 
300±i^ pc, using WISE photometry. 

For V1007 Her, we use the Wl band of WISE t o find 
a distance of 560±ig° PC. lGreiner. Schwarz fc Wenzell (|l998l ) 
give a lower limit of d > 250 pc. 

We obtain a distance of 5001^ pc for CW Hyi, again 
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from the Wl-band magnitude. Our estimate is slightly 
large compared to t he distance of ~ 250 pc reported by 
ISchwope etail l|2002h : however, the donor star features are 
not prominent in the optical spectrum, implying that the es- 
timate of the donor contrib ution to the optical flux is quite 
rough jSchwope et alj|2002h . 

Our distance estimate of 3801^0 P c f° r 
J1610. 1+0352 (also called RBS 1563) from WISE pho- 
tometry i s in g ood agreement the estimate of 320 pc that 
lAk et all (|2008T) obtain from 2MASS data. 

We use the 2MASS #-band magnitude of 1RXS 
J231603.9- 052713 (RBS 1973) to find a distance of 
46 °ii20 PC. lRodrigues et all l|2006l ) have estimated ~ 410 pc, 
in agreement with our value, from a model of the optical and 
infrared flux that includes cyclotron emission and a heated 
secondary. 

There is no published information on the distances of 
EQ Cet, EG Lyn, and RX J0859.1+0537 (RBS 734). We 
use WISE detections of these 3 systems in our distance es- 
timates. 



4-2.3 Systems for which we obtain only very weak 
distance constraints 

We are now left with 4 systems for which estimating dis- 
tances presents more serious problems. These are CD Ind, 
IW Eri, CV Hyi, and FH UMa. For CD Ind and IW Eri, the 
method discussed above produces distance estimates that 
are inconsistent with other information, while no near- or 
mid-IR detections are available for CV Hyi and FH UMa. 
We discuss these systems in turn below. In Section ^. 3. II we 
will show that these four very poorly constrained distances 
do not have an important effect on our space density esti- 
mate. This is because we are able to place lower limits on 
the distances, and, in all 4 cases, the resulting upper limit on 
the contribution of these systems to the total space density 
is small. 

Apparent magnitudes in the 2MASS and WISE bands, 
together with our assumptions about typical donor flux con- 
tributions in the IR bands, wo uld imply a distance o f roughly 
110 pc for CD Ind. However. ISchwope et all (| 19971 ) find no 
sign of the donor in their optical spectra, and give a lower 
limit of of d > 250, assuming a donor spectral type of M5. 

For IW Eri, WISE magnitudes would yield a distance 
of only about 80 pc, which is inconsistent with its non- 
detection in 2MASS. We measure K > 15.5 (from the 
faintest detections in the area around IW Eri), implying 
d > 150 pc. 

Both CD Ind and IW Eri have low- and high ph oto- 
metric states (|Schwope et al.lll997l ; ISchwope et alj|2002h . It 
is possible that CD Ind was observed in a high state by 
2MASS as well as WISE, while IW Eri was caught in a 
high state by WISE and in a low state by 2MASS. The 
high-state optical spectrum of CD Ind shows rising flux to- 
wards the red end, interpreted as cycl otron emission cause d 
by a weak magnetic field (B ~ 11 MGi lSchwope et al.|[l997l ). 
Therefore, CD Ind likely has unusually large cyclotron emis- 
sion in the IR, because of its unusually low magnetic field 
strength. 

For CV Hyi and FH UMa, we find lower limits of 
d > 300 pc and d > 320 pc, respectively, from WISE non- 



detections. iBurwitz et all (|l997n also report d > 300 pc for 
CV Hyi, based on its optical spectrum. 

We can derive conservative upper limits on the distances 
of these 4 systems from what would be implausible optical 
luminosities, or even position in the galaxy. However, since 
such limits are very weak, they are not useful here. It does 
seem that IW Eri is not ve ry distant, since it has a prope r 
motion of about 45 mas/yr i|Monet et al.ll2003l ; |Peterdl2008r ). 
iPattersonl (|201ll ) finds a tangential velocity of 39 ± 5km/s 
for a sample of normal non-magnetic CVs. If we assume that 
IW Eri comes from a population with the same age, we can 
tentatively estimate a distance of ~ 180 pc. 

For these 4 systems we will simply assume distance dis- 
tributions that are Gaussian in log(d), with 0"; O9 (d) chosen 
so that the 16th percentile corresponds to the lower limits 
given above. Of course this is not correct, but our space 
density estimate is not strongly affected (see Section l6.3.1l) . 

4.3 Possible bias in distance estimates 

Distance estimates used here may suffer from the well- 
known Malmquist l|Malmquisj 1 1924 in the case of pho- 
tometric parallax) and Lutz-Kelker |Lutz fc Ke lkcr 1973; 
for trigonometric parallax) biases. We are not concerned 
about Lutz-Kelker bias in the parallax measurements of 
EX Hya, TV Col and EF Eri, beca use the parallax errors 
are very small in the fi rst two cases (|McArthur et al.ll200ll ; 
iBeuermann et all 2003f). a nd because the bias was consid- 
ered bv lThorstensenl ( 2003) in the case of EF Eri. This leaves 
the possibility of Malmquist bias in the photometric paral- 
lax distances presented in Section [4.1,21 and 14.2.21 We have 
checked how large this effect cou l d be, in the same way as 
described in iPretorius fc Kniggel (|2012T ). We find that the 
bias, if present at all (se e the more detailed discussion in 
IPretorius fc Kniggdl2012T ). is in each case insignificant. 



5 X-RAY SPECTRA AND LUMINOSITIES 

Polars are luminous soft X-ray sourc es, and a large 
number was detected by ROSAT (e.g. H euermann et al.l 
1 19991 ; iThomas et al.l [199a ). Their X-ray spectra are usu- 
ally modelled as blackbodie s , with kT ~ 20 to 40 eV (e.g. 
Ramsay. Cropper fc Masonl Il996l; IThomas fc B cucrmann 



19981 ; Table 6.5 of IWarnerlll995h . Although a harder com- 



ponent is often present, in most polars it only contributes 
a significan t part of the flux in the RO SAT band during 
low states (Bcucrmann fc Thomas! Il993h . where even very 
nearby polars are too faint to be included in the RBS (see 
also Section |6.2.2|) . A handful of confirmed and candidate 
polars that show no soft X -ray component, even in the high 
state, is now known (e.g. iRamsav fc Cropped |2004| . 120071 ; 
iRamsav et alj|2009l ; IVogel et al.ll2008T l. However, these sys- 
tems do not appear to be intrinsically very common. Only 2 
are within 200 pc, and e.g. XMM-Newton and INTEGRAL 
and have not discovered many new polars (despite the near- 
est of the known hard polars, namely BY Cam, V1432 Aql, 
and V2301 Oph, being detected in the INTEGRAL/IBIS 
survey) . 

IPs have hard X-ray spectra, often modelled as intrin- 
sicall y absorbed thermal bremsstrahlung with kT ~ 30keV 
(e.g. lPatterson|[i994l ; iFrank. King fc Ramd["l985l ). In a few 
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IPs, a soft component is present (e.g. IHaberl fc Motchll 19951 ; 
Ide Martino et al1l2004l ). but our sample does not contain any 
of these "soft IPs". Depending on the accretion geometry, 
IPs can have very high intrinsic absorption (the intrinsic 
absorption also varies in individual systems). We take the 
intrinsic N H as 2 x 10 20 cm" 2 fo r EX Hya and DO Dra 
l|Richmanlll996l ; iMukai et al.ll2003l ). an d 4 x 10 22 cm" 2 for 
CC Scl, AO Psc TV Col, and EI UMa JPietsch et al II 19871 ; 
iRana et all 120041 ; iRamsav et all l2008l : IWoudt et al.l l2012h . 
with a covering fraction of 0.5 in all cases. 

We assume a kT — 30 eV blackbody spectrum for 
the polars, and a kT — 30keV thermal bremsstrahlung 
spectrum for the IPs in our sample. These single 
temperature, single component spectra are not very 
physical (e.g. iNorton fc Watsonl Il989l ; IMukai et all 12003 ; 
iBeuermann. Burwitz fc Reinschll2012l ). and they are not ex- 
pected to be good approximations for all sources in the sam- 
ple. The assumed intrinsic absorption for the IPs in the sam- 
ple are also rough approximations. Therefore, the values we 
give for Lx should be treated with caution. However, the 
space density estimate is much less sensitive to both the as- 
sumed X-ray spectra and the amount of intrinsic absorption 
we adopt for the IPs than are the estimated Lx values (we 
will return to this in Section \6. 3. 4[) . 

Since this is a high Galactic Latitude sample, in- 
terstellar absorption is low for all our sources. A few 
systems are at sufficiently large distances tha t the total 
Galactic Nh (as given by iKalberla et all [2005) is a good 
approximation. Many o f the polars have Nh estimates from 
X-ray s pectral fits (e.g . ISchwope et al.l 1997; Rcinsch et al 



iTinnev et all l|l993l ): 




Schwarz et all 120091; iBarrett. Singh fc Mitchel 



Ramsav, Cropper & Mason 
Pandel fc Cordoval I2005T 



19961; 



d. 



Schwope et al 



Ramsav et al 
2007f: 



Beuermann. Thomas fc PietschI Il99ll; IBurwitz et "al 



ISingh et~ 



19981 ; IBurwitz et all [l99 



iGreiner. Schwarz fc Wenzel 
For several systems, there 
are absorption measurements based on UV data (e.g . 
lArauio-Betancor et al.ll2005l; lLa Douslll99ll ; IVerbuntll 19871 ) . 
and for a few more, iBruch fc Engelllll994h give Ay estimates 
(we use the dust to gas ratio of iPredehl fc SchmittJ Il995l 
to obtain Nh)- In cases where we have no observational 
estimate of interstellar abso rption, we use the m odels of 
lAmores fc Lepind (|2005l ) and lDrimmel et all (|2003h . 

We assume 50% errors in estimates of interstellar col- 
umn densities, and give unabsorbed X-ray luminosities in 
the 0.12-2.48 keV band in Table □ 



[2 - (x) + 2xj + 2) 



The index j represents mCVs in the sample; Q is the solid 
angle observed in the survejQ, and Xj — dj\ saib\/h, with dj 
the maximum distance at which system j could have been 
detected (a function of its luminosity, the flux limit, and ex- 
tinction along all lines of sight covered by the survey) . This 
relation assumes that p does not vary with radial position in 
the Galaxy, and that the vertical density profile is exponen- 
tial. We assume scale-heights of h = 260 pc for short-period 
systems, and h = 120 pc for long-period systemfl To find 
dj, we need Nh as a function of distance and b. For this 
we assume an exponential vertical density profile for gas, 
with a scale-height of 140 pc, and normalized to give Nh 
corresponding to Ay — 0.8 mag/kpc in the Galactic Plan^D- 

The sum of the contributions of the 30 systems in the 
sample then gives the mid-Plane value of p (po = V^j)- 
For the moment ignoring all errors, and simply assuming 
that the best-estimate distance, count rate, and Nh are the 
true values, we list the fractional contribution that each sys- 
tem in our sample makes to the total space density {pj/po) 
in the 8th column of Table [1] This shows that our p mea- 
surement is not dominated by 1 or 2 objects. 

In order to correctly sample the full parameter space 
allowed by the data, and hence find the error on po, we 
compute its probability distribution function using a Monte 
Carlo simulation that calculates po for a large number of 
mock samples. The mock samples are created by drawing a 
distance, Nh , and count rate for every observed system from 
the appropriate distribution. These values are used to calcu- 
late each 1/Vj, which is then weighted by a factor fij, drawn 
from the probability distribution of the number of sources 
belonging to the population (corresponding to a particular 
observed system) that one expects to detect in the R BS (see 
IPretorius et~aT]|2007bl and IPretorius fc Kniggell2012l ). 



6.2 Results 

We have carried out this calculation for the whole observed 
sample of 30 mCVs, as well as for several sub-samples (po- 
lars, IPs, short- and long-period systems, and short-period 
polars and long-period IPs). Table [2] summarizes the results. 



6 CALCULATING THE SPACE DENSITY 
6.1 The method 

We describe the calculation that gives the space den- 
sity, together with the uncertainty on that estimate 
only briefly. More detaile d di scussions may be foun d in 
IPretorius fc Kniggel |2012l ) and IPretorius et all (|2007bl ). 

The effective observed volume of the RBS depends on 
the spacial distribution of CVs in the Galaxy, and on the 
X-ray luminosities of systems in the sample (since the sur- 
vey is flux limited, rather than vo lume limited) . It is found 
using the relation given by e.g. IStobie et al] (|l989l ) and 



5 Because b is variable over f2, we compute each Vj as a sum 
over smaller solid angles, SQ. The terms in this sum are slices 
subtended by 2° in b. These 5Q are sufficiently small that the 
error introduced by b varying over 5Q is negligible. 

6 The reason for this is that short- and long-period CVs are ex- 
pected to be old and young populations, respectively. It is only 
a crude approximation, since systems can form at short period, 
and since polars at all periods are believed to evolve slowly. Our 
results do not change significantly if we simply assume the same 
h for all systems (see Section l6.3.2D 

7 Since the survey covers a large area, at high galactic latitudes 
(all |6| > 30°), and the maximum distances we deal with are large 
compared to the typical size of inhomogencitics in the ISM, this 
smooth model for the density of gas is justified (for the most part, 
but see Section l6. 2.21 1. We explore the sensitivity of our p estimate 
to the assumed mid-Plane gas density in Section 16.3.31 
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Table 2. Space density estimates for all mCVs, IPs and polars 
separately, short- and long-period mCVs, and short-period polars 
and long-period IPs. UW Pic and J0859 + 05 are excluded from 
the samples split by P or i,, because their orbital periods fall in the 
period gap. In addition to the values inferred from the observed 
RBS sample, we list the estimates that would be obtained if polars 
are undercounted by a factor of 2, because of a high-state duty 
cycle of only 0.5. 



1.5xl0 6 



Sample 



Po/(10~ 7 
No assumed low states 



pc 3 ) 
Polar duty cycle 0.5 



All mCVs 


o 9 +3.9 


13 +6.3 
10 -3.8 


Polars 


4 9+2-7 


gg+5.4 
a -°-3.1 


IPs 


9 7 +2.4 
Z -' -1.2 




short-P oro mCVs 


4.6 +3 / 2 


o o+5.7 
°' d -3.0 


long-P,,,.;, mCVs 


2.8t\i 


3 R+2-2 
0.0_ 1 4 


short-P or f, polars 


o c+2.5 
a - D -1.9 


~ 9+5.O 
' - Z -2.7 


long-P ort IPs 


i 7 +i.6 
i -'-0.9 





6.2.1 The probability distribution function of po 

The distribution of po values, normalized to give a proba- 
bility distribution function, from the simulation including 
all magnetic systems is shown in Fig. [2] The mode, me- 
dian, and mean of the distribution are marked by solid lines 
at 7.1 x 1CT 7 , 8.2 x 1CT 7 , and 9.2 x 10~ 7 pc~ 3 , while the 



dashed lines at 5.8 x 10 and 1.2 x 10 pc show a l-o 
confidence interval (the 16th and 84th percentile points of 
the distribution). In other words, our estimate for the mid- 
plane space density of mCVs is 8^2 x 10 _7 pc -3 . We find 
that the large errors in many of the distances dominate the 
total uncertainty in po, although the small sample size and 
Nh errors also contribute significantly. 

The inset in Fig. [5] shows the probability distribution 
functions of po for the whole mCV sample again, as well 
as for IPs (dotted histogram with lowest mode) and polars 
(dashed histogram) separately. We find space densities of 



3± 2 x 10~ 7 pc~ 3 for IPs and 51?, x 10~ 7 pc~ 3 for polars. 



6. 2. 2 Upper limits on the space density of an undetected 
population 

The estimate above assumes that the detected population of 
mCVs is representative of the true underlying population, 
in the sense that it contains at least 1 of the faintest IPs and 
polars that occur in the intrinsic population (but not that 
it contains faint and bright systems in proporti on to their 
intrinsic incidence; see IPretorius fc KniggellioTj ). Since the 
effective volume of the survey is smaller for fainter Lx , it is 
possible that even a large population of sources at the faint 
end of the luminosity function can go completely undetected. 
Here we place limits on the sizes of faint populations of 
polars and IPs that could escape detection in the RBS. 

We again perform a Monte Carlo simulation with the 
same simple Galaxy model for stars and gas as described in 
Section [671] above (but, for simplicity, now assuming a single 
scale height of 260 pc for all mCVs). A model population of 
mCVs, all with the same Lx, and with a spectrum appro- 
priate to either polars or IP£j, is distributed in the model 

8 For polars, we again assume a kT = 30 eV blackbody spectrum, 
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Figure 2. The po distribution for all mCVs, resulting from our 
simulation. Solid lines mark the mode, median, and mean at 7.1 X 



Dashed lines show a 1- 
-6 „„-3 



10" 7 , 8.2 x 10~ 7 , and 9.2 x 10~ 7 pe- 
er interval from 5.8 X 10~ 7 to 1.2 x 10~ t, pc~ ;i . The probability 
distribution functions shown in the inset are (with modes from 
high to low po) for the whole mCV sample, polars alone (dashed), 
and IPs alone (dotted). In the inset, solid lines at 8.2 x 10 -7 , 
4.9 X 10 -7 , and 2.7 x 10~ 7 pc -3 mark the medians of the three 
distributions. 



Galaxy, in order to find the value of po for which the pre- 
dicted number of detected systems is 3 (so that detecting 
such systems is a 2-er result). We do this for a range of Lx 
below the faintest values found for the observed sample. 

Fig. [3] shows the maximum allowed po as a function of 
Lx, separately for possible undetected polar and IP popu- 
lations. The limits from the simulation are plotted as bold 
histograms, and the fine curves are fits to the data, given by 



Pmax = 1.02 x 10- 5 (Ljf/10 30 ergs- 
for IPs, and 

Pmax = 4.01 X III ":/.X II) .'I- S 



pc 



pc 



for polars. p m ax is only the upper limit on the mid-plane 
space density of an undetected population, and does not 
include the contribution from the observed systems. Thus, 
as a specific example, a population of undetected polars with 
a space density as high as 5x the measured p po i ar must have 



po = 5 x pip if it consists of systems with X-ray luminosities 
fainter than 5 x 10 30 ergs _1 . 

Because the two relations given above involve Lx , which 
depends quite sensitively on our assumed X-ray spectrfl 
they are of only limited use (this is discussed further in Sec- 
tion [OTj}. 

This is also an instance where our smooth model of 
the ISM is not a valid approximation. This is because, e.g., 

and for IPs, we assume a kT = 30keV thermal bremsstrahlung 
spectrum with a partial covering absorber with covering fraction 
of 0.5 and N H = 2 x 10 20 cm" 2 . 

9 For example, adopting blackbody temperatures of 10 and 50 eV 
for a moderately absorbed polar (Nh = 2 x 10 2t) cm - 2 ), yield 
values of Lx that differ by a factor of 20. For the hard spectra 
of IPs, Lx is less sensitive to the assumed temperature of the 
bremsstrahlung spectrum, but does vary widely over an intrinsic 
Nh range as large as is observed. 
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polars with Lx = 10 ergs can only be detected out to 
~ 40 pc, which places them in side the 'Local Bubble' (e.g. 
iFrisch. Redfield fc Slavinll201lh . However, since our model 
probably gives too high absorption at small distances, it 
only means that, at least for polars, the upper limit on po 
is conservative at the faint Lx end. In the case of IPs, the 
structure of the local ISM matters less, since their spectra 
are less affected by interstellar absorption. 

Polars probably do not experience a large range in sec- 
ular M over the course of their evolution (because M is a 
relatively flat function of P or b when the only AML mech- 
anism is GR). This means that we do not expect the po- 
lar luminosity function to rise towards the faint Lx end, 
beyond the luminosities that we are sensitive to. However, 
most (possibly all) polars switch between low- and high M 
states on shorter timescales (e.g.. [Ramsay et al.ll2004bT ). In 
the low photometric state, even very nearby polars are too 
faint to be included in the R BS. An example is AR UMa. It 
is at a distance of only 86 pc (|Thorstensen. Lepine k, Sharal 
20081'). and is a bright soft X-ray source in its high state (e.g. 
Remillard et al .11 19941 ). but was so faint during the RASS 
that it is not even included in the faint source catalogue. 
The deep low states of polars imply that in a single-epoch 
X-ray survey (which is what the RASS was over most of 
the sky) only a fraction of the local polar population that is 
equal to the high-state duty cycle can be detected. In other 
words, our assumption that the polar sample detected in 
the RBS is representative of the intrinsic population, is not 
strictly valid, since it does not include low-state systems. 
I Ramsay et al.l (|2004bh find that polars spend roughly half 
their time in low states, implying that our estimate of p po iar 
is probably a factor of 2 too lovQ In the final column of 
Table [21 we list the space density estimates that we obtain 
when allowing for a high-state duty cycle of 0.5 for polars 
and assuming that polars are undetectable their low states. 



6.3 Sensitivity of the results to uncertain 
assumptions 

Here we examine the sensitivity of our space density esti- 
mates to some of the more uncertain assumptions we made 
in calculating it. We will consider the possible impact of 
the systems with the most poorly constrained distances, the 
scale-heights we adopt for different CV populations, the as- 
sumed amount of interstellar absorption, the completeness 
of the sample at lower |6|, and the assumed X-ray spectra. 
By investigating an extreme range of assumptions, we show 
that possible systematic errors affect these space density es- 
timates to less than roughly a factor of 2. 

6.3.1 The most uncertain distance estimates 

Four members of the mCV sample, CV Hyi, FH UMa, IW 
Eri, and CD Ind, have very poorly constrained distances (see 
Section [4.2.3|> . The lower limits on distances are, however, 



Note that the high-state d uty cycle of polars is not 
very w ell constrained; see e . g. iHess man, Gansicke. fc Matteil 
<2000h. lArauio-Betancor et al.l J2005t >. IWu fc Kisj j200Sf) , and 
iBreedt et al.l J2012T) . 
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Figure 3. The upper limit on the mid-plane space density as 
a function of X-ray luminosity for an undetected population of 
polars and IPs. The data from the simulations are shown as a 
bold histogram, and fits are over-plotted as a finer lines. Note 
that the assumed X-ray spectra of polars and IPs are different, 
hence the different slopes. 

always robust, and allow us to determine that these sys- 
tems make at most a small contribution to the overall space 
density. Specifically, if we assume in each case the smallest 
allowed distance, the median of the po distribution would 
be at 8.7 x 10~ 7 pc -3 , while neglecting the contributions of 
these 4 systems yield 6.6 x 10 _7 pc -3 . This is within the er- 
rors of our best estimate of 8jl 2 x 10 -7 pc -3 . At its minimum 
allowed distance, IW Eri would account for 0.1 of the total 
space density, while the other 3 would all contribute less. 

6.3.2 The scale height of mCVs 

We have assumed scale-heights of 260 and 120 pc for short- 
and long-period systems, respectively, as a crude approxi- 
mation to the expectation that short- and long-period CVs 
belong to populations of different typical ages. If we instead 
set h to 260 or 120 pc for all systems, we obtain po distribu- 
tions with medians at 6.3 x 10 _7 pc~ 3 and 14 x 10~ 7 pc~ 3 , 
respectively. Therefore, although a plausible range of scale- 
heights shift the po distribution considerably, one would ob- 
tain estimates consistent with our result for other reasonable 
assumptions regarding the Galactic distribution of mCVs. 

We have also checked that the ^-distribution of the RBS 
mCV sample is consistent with the Galaxy model we use. 
In this calculation, we use the same model for the density 
of stars and gas, and assume the same X-ray spectra as be- 
forcE 

After imposing the RBS flux- and |6|-limits, we use 
a Kolmogorov-Smirnov (KS) test to compare model and ob- 
served ^-distribution. We find that the whole sample of 30 
observed systems, as well as subsamples consisting of short- 

11 Since the RBS is a high Galactic latitude, flux-limited survey, 
the observed sample does not have the same z-distribution as the 
underlying population. Generating a model z-distribution there- 
fore involves an assumption regarding the luminosity function of 
the intrinsic polar and IP populations. We have assumed intrin- 
sic distributions that are Gaussian in log(Ljf ), and experimented 
with several values for the average and standard deviation. 
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and long-period systems, have ^-distributions that are con- 
sistent with both h — 260 and 120 pc. In other words, our 
sample is too small, and the distance errors are too large, to 
distinguish between these two choices of scale height. 



6.3.3 Interstellar absorption and completeness of the 
sample 

Our treatment of interstellar absorption is very simple, and 
could be a concern, considering the soft, easily absorbed 
spectra of polars. As noted before, the unrealistically smooth 
model of the interstellar medium is not an important short- 
coming (since we are using it to find a survey volume which 
covers a large fraction of the whole sky, at relatively high 
|6|, and since even the smaller maximum distances involved 
in the calculation are large enough to implying that we can 
average over regions of low and high absorption) . The total 
amount of absorption, on the other hand, could have a large 
effect on the space density estimate, because higher absorp- 
tion (along all lines of sight) reduces the survey volume. 

We have set the mid-plane density of gas to produce 
Ay — 0.8 mag/kpc for b — 0°. A wid e range in average mid - 
plane extinction has been reported. iDrimmel et alj (|2003h 
give Av — 0.7 mag/kpc for lines of sight near the Galactic 
Plane in the inner disc of the Galaxy, and Ay ~ 0.5 mag/kpc 



in the outer Galactic disc. The model of lAmores fc Lepinel 
(|2005h produces Av ~ 1 mag/kpc, and lVergelv et all (|l998T ) 
find 1.2 mag/kpc in the Galactic Plane. Repeating our space 
density calculation with the density of interstellar gas nor- 
malized to give values of 0.5 to 1.2 mag/kpc in the Plane, 
yields po distributions with median values of5.8xl0~ 7 and 
12 x 10 _7 pc -3 , respectively^, This is consistent with our 
best estimate of 81^ x 10~ 7 pc~ 3 . 

We can also check for evidence that the completeness 
of our sample decreases at lower \b\, as one might expect 
from increasing absorption along lines of sight closer to the 
Galactic Plane. To do this we simply create subsamples with 
different \b\ cutoffs. Compared to p = 8l* x 10 -7 pc -3 from 
the whole sample of 30 systems (|6| > 30°), we find median 
po values of 8.1 x 10~ 7 pc~ 3 and 5.6 x 10~ 7 pc~ 3 for \b\ > 40° 
(17 systems) and |6| > 50° (10 systems), respectively. These 
number are consistent to within the errors, and if there is a 
trend, it goes in the opposite direction than expected if the 
sample is less complete at lower \b\. 



6.3.4 The assumed X-ray spectra 

A final source of systematic uncertainty associated with the 
space density estimate is our assumed X-ray spectra. We 
have used simple, single component spectra for all systems 
(bremsstrahlung spectrum for IPs, and blackbodies polars). 
Here we will check the possible impact of a range of X-ray 
spectra on the space density estimates. 

The space density calculation uses the maximum dis- 
tance at which a given mCV in the observed sample could 
have been detected, since this determines the survey volume 
(for a population of systems with the same Lx). For a given 



12 pip on its own is much less sensitive to this, because of the 
harder spectra of IPs, and because intrinsic absorption reduces 
the importance of interstellar absorption. 
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Figure 4. The ratio of the values of 1/Vj obtained by assuming 
kT = 10 and 50 eV blackbody spectra, as a function of the count 
rate at which the system is detected. The 4 curves are for different 
true distances. The count rates of 4 of the polars in our sample 
are indicated on the curves corresponding to their best-estimate 
distances. The largest effect of different spectra on pj occurs for 
a system that lies at a small distance compared to the maximum 
distance it could have been detected at (since this implies a large 
amount of absorption between d and d ma x). The most extreme 
example in the RBS is EF Eri, where its contribution to the total 
space density would be a factor of about 2.2 larger for an assumed 
kT = 10 eV than for an assumed kT = 50 eV spectrum. 



object, the maximum distance depends on the ratio of Fx 
to the flux limit, which may as well be expressed as a ratio 
of observed count rate to limiting count rate. This means 
that the assumed X-ray spectral shape has relatively little 
influence on the space density calculation (only the inter- 
stellar absorption depends on the X-ray spectrum, and this 
only matters to the extent that Nh differs between the true 
and maximum distance, and that the slope of absorption as 
a function of Nh differs for different spectra). 

To show this more explicitly, if y is the ratio of count 
rate (c) to X-ray flux (i.e. y — c/Fx, a function of the 
spectrum, the instrument response, and Nh), then the dis- 
tance is d = y/yLx/^Trc, so dmax = \/y m axc/yci im , where 
Umax = y(dmax) and ci im is the limiting count rate. Then, for 
d ~ dmax (the case of a system detected at close to the lim- 
iting count rate), and for large d (Nh at d already at about 
the total value for the Galaxy), y max ~ y, implying that 
the shape of y(Nn) does not matter, and different spectra 
give the same dmax (and therefore 1/Vmax)- Otherwise, the 
difference in dmax resulting from assuming different spectra 
depends on the difference in slope of y(Nn) for those spectra. 
We illustrate this in Fig. [4] where we use blackbody spec- 
tra with kT — 10 and 50 eV. The figure shows that, while 
different assumed spectra imply very different values of Lx , 
the space density contribution of a given system changes by 
at most a factor of ~ 2, with a much smaller effect in most 
cases. The effect is also smaller for spectra appropriate to 
IPs than for the soft spectra of polars. 

For the whole sample of polars, we would obtain median 
po values of 8.1 x 10~ 7 pc -3 and 4.5 x 10~ 7 pc~ 3 assuming a 
kT — 10 and 50 eV blackbody spectrum, respectively, com- 
pared to Ppoiar = 5^2 x 10~ 7 pc -3 for our favoured 30 eV 
spectrum. Note that while a kT = 30 eV blackbody spec- 
trum is unlikely to be a good approximation for all the polars 
in our sample, we also do not expect all to be better rep- 
resented by either much lower or much higher temperature 
blackbodies. 
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pip is less sensitive to the assumed spectrum, because 
the difference in absorption between different hard, intrin- 
sically absorbed spectra is small, even for a large range in 
assumed temperature. Our best model is a kT = 30keV 
thermal bremsstrahlung spectrum, with intrinsic column 
densities of 2 x 10 20 cm" 2 for EX Hya and DO Dra, and 
4 x 10 22 cm" 2 for CC Scl, AO Psc, TV Col, and EI UMa, all 
with a covering fraction of 0.5. We will separately consider 
different bremsstrahlung temperatures and intrinsic absorp- 
tion. With kT = 10 and 50 keV, respectively, we obtain 
median po values of 2.8 x 10~ 7 pc~ 3 and 2.6 x 10~ 7 pc~ 3 
for the sample of 6 IPs, compared to our best estimate of 
3± 2 x 10" 7 pc" 3 . Assuming kT = 30keV for IPs, but with 
either high (Nh — 4 x 10 22 cm -2 and a covering fraction of 1 
for all systems) or no intrinsic absorption, we find po distri- 
butions with medians at 2.5 x 10 -7 pc -3 and 2.8 x 10~ 7 pc -3 , 
respectively. 

The upper limits on the space densities of hypotheti- 
cal hidden populations of IPs and polars are of course also 
affected by the assumed X-ray spectrum. Since we find the 
upper limits as a function of the X-ray luminosity of the 
undetected population (and Lx, for the detected as well 
as "hidden" systems, depends on the assumed X-ray spec- 
trum), we present the effect of assuming different X-ray spec- 
tra as the ratio of p m ax and the best-estimate measured 
po, as a function of the ratio of Lx of the hidden popula- 
tion and the faintest detected system. In other words, we 
fix Ppoiar — 5 x 10 _7 pc -3 and pip — 3 x 10 _7 pc -3 , and 
then find p m ax, P oiar and p ma x,ip, as well as L x of the hid- 
den population of IPs and polars and of our faintest ob- 
served IP and polar with different assumed spectra. The 
result is shown in Fig. [5] In the case of polars, a popula- 
tion with a space density as high as 5x the measured p po iar 
(i.e., p max ~ 4 x 10~ 6 pc~ 3 must have L x < 0.06x the 
luminosity of EQ Cet, the faintest polar in the RBS sam- 
ple. An undetected IP population must have Lx J$ 0.2 x 
the luminosity of the faintest detected system (CC Scl) in 
order for its space density to be as high as 5 x pip (i.e., 
Pmax — 1.4 x 10 -6 pc -3 ). These ratios are not very sensitive 
to the spectrum we assume (see the size of the hatched areas 
in Fig. II}. 



7 DISCUSSION 

Having derived space density estimates for polars and IPs, 
let us revisit the three key questions that provided the main 
motivation for our study (see Sections [1] and [2} . We note 
from the outset that our goal here is merely to take a brief 
look at some of the implications of our results; a full analysis 
of these issues is well beyond the scope of the present paper. 



7.1 The relationship between intermediate polars 
and polars 

As noted in Section [2 a natural explanation for the lack of 
short- P or b IPs is that long-period IPs evolve into polars be- 
low the period gap. If one assumes that long-period IPs are 
the sole progenitors of short-period polars, and that all IPs 
synchronize once they have crossed the period gap, then the 
ratio of the space densities of long-P or (, IPs and short- P or b 




Figure 5. Pmax,Ip/pIP and p ma x,polar/ Ppoiar as a function of 
the ratios of the luminosity of a hidden population and the lumi- 
nosity of the faintest detected IP/polar (CC Scl/EQ Cet). The 
bold curves are for our best model (a 30 eV blackbody spectrum 
for polars and partially absorbed 30 keV thermal bremsstrahlung 
spectra for IPs) . The hatched areas show the sensitivity to the as- 
sumed spectra. For polars, this corresponds to the area between 
the curves obtained by assuming respectively kT = 10 and 50 eV 
blackbody spectra in calculating pmax, polar and X-ray luminosi- 
ties for both the undetected population and EQ Cet (while p po i ar 
is fixed at 5 X 10 -7 pc - 3 ). For IPs, the hatched area corresponds 
to kT = 10 and 50 keV bremsstrahlung spectra. 

polars (pip t i p and p po i ar ,sp) should simply reflect their rela- 
tive evolutionary time-scales. For the purpose of our qualita- 
tive discussion, let us throw caution to the wind and assume 
that the numbers in Table [5] are reasonable approximations 
to the total space densities, even though our samples are 
flux-limited and even though polars and IPs have quite dif- 
ferent X-ray luminosity functions. 

We then find that the observed logarithm of this ra- 
tio is log (ppoiar.sp/ pip,i P ) = 0.32 ± 0.36 (this becomes 
log {ppoiar.sp I ' Pi p,i p ) = 0.63 ± 0.36, if we assume a 0.5 high- 
state duty cycle for polars). The ratio itself is therefore ~ 2, 
but this is accurate to only about a factor of two. If the evo- 
lution of long-period IPs is really driven by MB, while that 
of short-period polars is driven solely by GR, the evolution- 
ary time-scale of the latter is expected to exceed that of th e 
former by at least a factor of >5 (e.g. iKnigge et al]|201ll '). 
This is larger than the ratio of the inferred space densities, 
but still completely consistent with it, given the rather large 
statistical errors. In fact, at 2-cr, the uncertainties are large 
enough to encompass both ratios exceeding 10 and ratios 
less than unity. This means that, with the currently avail- 
able space density estimates for polars and IPs, we cannot 
place strong constraints on the evolutionary relationship be- 
tween the two classes. Nevertheless, it is interesting to note 
that the simplest possible model, in which short-period po- 
lars derive from long-period IPs, is not ruled out by their 
observed space densities. 
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7.2 The intrinsic fraction of magnetic CVs 

We can also combine our measurement of the space density 
of m CVs (polars as well as IPs) with that of non-magn etic 
CVs jPretorius fc Kniggel 1201a iPretorius etU] l2007bh in 
order to obtain an estimate of the intrinsic fraction of mCVs 
amongst the Galactic CV population (f m cv)- This again in- 
volves throwing caution to the wind to some extent, since 
the samples involved are flux-limited, and the populations 
being compared have different X-ray luminosity functions. 
Thus it is quite possible that our estimate of the mCV frac- 
tion is still affected by selection biases. 

Keeping this caveat in mind, we find that Iog(/ m ov) = 
— 0. 80^0 gg, i.e fmcv — 16%, to within a factor of 2 (or 
log(/ m cv) — —0.63^0 33, if half of all polars are in the low 
state at a given time, and thus undetected in the RBS). 
This is consistent with the ra w incidence of mCVs in the 
known CV sample (~20%; see iRitter fc Kolbll2003h . How- 
ever, more importantly, it is also consistent, within our con- 
siderable uncertainties, with the fraction of isolated WDs 
that are strongly magnetic (~ 10%). In fact, it seems likely 
that our X-ray-selected CV sample is more complete for 
mCVs than it is for non-magnetic CVs. If so, our estimate 
of fmcv should be considered as an upper limit. Thus the 
incidence of magnetism is not obviously enhanced amongst 
CV primaries compared to isolated WDs. 



7.3 Galactic X-ray Source Populations 

As noted in Section [2] mCVs - and particularly IPs - have 
been suggested to be the dominant X-ray source popula- 
tions above Lx — 10 31 ergs _1 in a variety of Galactic envi- 
ronments. We can use our new space density estimates to 
check whether mCVs can plausibly account for the number 
of sources seen in surveys of these environments. 

Let us take the Galactic Centre regio n as a n example. 
The deep Chandra survey of lMuno et all ( 2009) covers an 
effective area of ~ 10 _3 deg 2 down to Lx — 10 31 ergs -1 and 
includes ~ 9000 sources. For our order of magnitude esti- 
mate here, we will ignore subtleties like the flux/luminosity- 
dependent survey area and simply ask if it is plausible that 
the majority of these sources may be IPs. 

Given that the stellar density distribution is highly 
peaked towards the Galactic Centre, let us approximate 
the volume covered by the survey as a sphere of radius 
R ~ 150 pc. The space density of X-ray sources in the Galac- 
tic Centre is then of order px,GC ~6x 10 -4 pc -3 , while the 
local space density of IPs is pip ~ 3 x 10 -7 pc -3 . However, 
the stellar space density in the Galactic centre is ~ 70 pc -3 , 
while it is only ~ 0.044 pc -3 in the solar neighborhood (e.g. 
iHong et alj2009h . We thus find that there is roughly 1 X-ray 
source per 100,000 stars in the Galactic centre, and roughly 
1 IP per 200,000 stars in the solar neighborhood. At the level 
of precision to which we are working here, these numbers are 
identical. We thus conclude that IPs remain a viable expla- 
nation for most of the X-ray sources seen in the Galactic 
Centre. 

It should be obvious that the calculation above is not to 
be taken too seriously. Its purpose is merely to illustrate how 
our measurement of the space density of mCVs, in general, 
and IPs, in particular, relates to recent X-ray surveys in 
a wide variety of Galactic environments. A correct analysis 



would have to account in much more detail for the properties 
of the various surveys. This is worth doing, but beyond the 
scope of the present paper. 



7.4 Outlook 

The rough calculations above show that open questions con- 
cerning the evolution and Galactic abundance of mCVs can- 
not yet be conclusively answered. This highlights a funda- 
mental problem: given the limited size of the existing flux- 
limited CV samples, and the low precision of most avail- 
able distance estimates, it is currently impossible to measure 
space densities to an accuracy much better than a factor of 
~ 2. To make matters worse, we often need ratios of space 
densities for various sub-populations in order to test evo- 
lutionary models, which necessarily suffer from even lower 
precision. 

This situation should improve dramatically over the 
coming years. Surveys with eROSITA will reach flux lim- 
its around 2 orders of magnitude deeper than the RBS and 
will y ield large X-ray-selected mCV samples (e.g. ISchwopel 
|2012| ). while Gaia will provide accurate distance measure- 
ments for a large number of CVs. 



8 CONCLUSIONS 

We have used a complete, purely X-ray flux-limited sample 
of 30 mCVs from the RBS to place constraints on the local 
space density of these systems. Our conclusions are listed 
below. 

(i) Assuming that the sample used here is representative 
of the intrinsic population (in the sense that the RBS de- 
tected at least one IP and one polar at the faintest ends 
of the luminosity functions of those populations), we obtain 
a mid-plane space density of Si 4 , x 10~ 7 pc -3 for mCVs. 
For the two distinct types of mCVs, we find p po iar = 
5±l x 10" 7 pc" 3 and p IP = Zt\ x 10~ 7 pc~ 3 . If we as- 
sume that polars are detectable in X-rays for only 50% 
of the time, we obtain p vo i ar = l^g 3 x 10 -6 pc -3 and 
p m cv = l-3±S:e x 10" 6 pc- 3 . 

(ii) We have calculated the maximum sizes of hypothet- 
ical faint populations of IPs and polars that are consistent 
with their non-detection in the RBS, as a function of the 
X-ray luminosity of the undetected populations. If an unde- 
tected population of polars with a space density 5x as high 
as the space density we infer from detected systems exists, 
then those systems must have Lx J$ 0.06 x the luminosity 
of the faintest polar in the RBS sample (EQ Cet). In the 
case of IPs, an undetected population with Lx J$ 0.2 x the 
luminosity of the faintest detected system (CC Scl) can have 
a space density as high as 5 x the value we measure from the 
detected IPs. 

(iii) The ratio of the space density of short-period polars 
to long-period IPs is 2+ 3 (or 4^, assuming a 50% high- 
state duty cycle for polars). Within the large errors, this is 
consistent with the very simple hypothesis that (the major- 
ity of) long-period IPs evolve into short-period polars, and 
that this accounts for the whole population of short-period 
polars. 



14 M.L. Pretorius, C. Knigge and A.D. Schwope 



(iv) Our estimate of the intrinsic fraction of mCVs is con- 
sistent with the fraction of magnetic systems in the known 
CV sample (~ 0.2). However, with existing data, this frac- 
tion cannot be measured to high enough precision to rule 
out an incidence of mCVs as low as ~10%. It is therefore 
not clear whether the fraction of strongly magnetic WDs is 
higher in CVs than in the single WD population. 

(v) When the local space density of IPs is scaled to the 
density of stars in the Galactic Centre, it is sufficiently high 
to account for the number of bright (Lx > 10 31 ergs -1 ) X- 
ray sources detected in that region. 
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